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ABSTRACT 



(N 

Aims. Current mass-loss rate estimates imply that main sequence line-driven winds are not sufficient to strip away the H-rich envelope 
^0 ' to yield Wolf-Rayet (WR) stars. The rich transitional population of the young massive cluster Westerlund 1 (Wd 1) provides an ideal 

laboratory to observe and constrain mass-loss processes throughout the transitional phase of stellar evolution. 

Methods. We present an analysis of deep radio continuum observations of Wd 1 obtained with the Australia Telescope Compact 
Array at four frequency bands that permit investigation of the intrinsic characteristics of the radio emission. 
^ . Results. We detect 18 cluster members, a sample dominated by the cool hypergiants, with additional detections amongst the hotter 

( ' OB supergiants and WR stars. The radio properties of the sample are diverse, with thermal, non-thermal and composite thermal/non- 

' thermal sources present. Mass-loss rates determined for stars with partially optically thick stellar winds are ~ 10"^ Mq yr"' across 

, all spectral types, insufficient to enable the formation of WRs during a massive star lifetime, and the stars must undergo a period of 

greatly enhanced mass loss. The sgB[e] star W9, the brightest radio source in Wd 1, may provide an example, with a current mass-loss 
rate an order of magnitude higher than the other cluster members, and an extended nebula interpreted as a wind from an earlier epoch 
with a density ~ 3x the current wind. Such an envelope structure in W9 is reminiscent of luminous blue variables, and one that shows 
^ ' evidence of two eras of high, possibly eruptive mass loss. Surprisingly, three of the OB supergiants are detected, implying unusually 

dense winds, though they are embedded in more extended emission regions that may influence the derived parameters. They also may 
have composite spectra, suggesting binarity, which can lead to a higher flux than expected from a stellar wind. Spatially resolved 
nebulae are associated with three of the four RSGs and three of the six YHGs in the cluster, which are due to quiescent mass loss 
\l I rather than outbursts. The extended nebulae of W20 and W26 have a cometary morphology, implying significant interaction with 

, either the intracluster medium or cluster wind. For some of the cool star winds, the ionizing source may be a companion star though 

the cluster radiation density is sufficiently high to provide the necessary ionizing radiation. Five WR stars are detected with composite 
' spectra, interpreted as arising in colliding-wind binaries. 
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1 . Introduction than just stellar evolution. For example massive stars are thought 

to predominantly form in stellar aggregates, where they drive 

Recent work strongly suggests that canonical mass-loss rates for cluster winds, which are a major source of mechanical feed- 

O stars need to be revised downwards to accommodate the effect ^^^k and chemically processed material into the wider galactic 

of wind clumping (FuUerton et al. 2006; Mokiem et al. 2007). environment, in turn driving star formation and galactic evolu- 

Such main sequence (MS) mass-loss rates are insufficient to re- ^^^^ indeed the feedback from large populations of Super Star 

move the H-rich mantle of the star prior to it becoming a Wolf clusters (SSC; Mtotai ~ 10^ - lO^M^) may drive galactic scale 

Rayet (WR), shifting the burden of mass loss onto the short lived outflows, which, if present in Dwarf galaxies, may be sufficient 

transitional phase of stellar evolution. This phase is populated by ^o strip them of their interstellar medium, preventing subsequent 

a wide variety of highly luminous, hot supergiant B[e] and lu- generations of star formation (e.g. Westmoquette et al. 2007). 

minous blue variable (LBV) stars, and cool Yellow Hypergiant ^ u i » ..u i • \ir » i j i ^wr » i j 

. . • /r^o^N . . One such cluster in the Galaxy IS Westerlund 1 (Westerlund 

(YHG) and red supergiant (RSG) stars. However the exact path Tr,/zi\ u f» wrj i e u- u u » j 

^_ , , r ■ r - ■ 1961), hereafter Wd 1, tor which photometric and spectroscopic 

or an O star through this transitional zoo as a function of ini- , , . . , fu i ju . 

°, , , ,. , observations suggested a unique population of both cool and hot 

tial mass is currently poorly understood, while the short-lived /r> » i i mn » i j Tnon\ n »i 

-'. ^ ^ , ' . ^ supergiants (Borgman et al. 1970; Westerlund 1987). Recently, 

epochs, and hence rarity, of such stars complicates efforts to con- . •, ■ -.• i j ir. u u j »u 

^ . ' . , , ^ ,,.r. detailed optical and near-IR observations have connrmed these 

strain their properties such as mass-loss rate and lifetime. , , ,,,, , ^ , » 

^ ^ results and revealed Wd 1 to be even more extreme than pre- 

A better understanding of such short-lived phases in the life yjo^giy anticipated, containing a large population of post-MS 
cycle of massive stars is crucial to areas of astrophysics other gj^rs with representative members of all evolutionary stages: OB 
supergiants and hypergiants, RSGs, YHGs and WRs (Clark & 
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respectively). Indeed, Wd 1 contains 6 YHGs, more than 50% 
of the currently known population in the Galaxy, as well as 
one of the largest WR populations of any cluster in the Galaxy 
(Crowther et al. 2006). With a cluster mass of -IO^Mq (COS), 
Wd 1 is directly comparable to the SSCs observed in external 
galaxies such as M82 and thus represents a relatively nearby ex- 
ample that provides a valuable opportunity to study the proper- 
ties, evolution and interaction of massive stars in their 'natural' 
environment. 

Radio continuum observations are a long established tool for 
estimating mass-loss rates for early-type stars. As part of a pro- 
gramme to accomplish this for classical Be and B[e] stars, (Clark 
et al. 1998, hereafter C98) imaged Wd 1 at radio wavelengths 
and found two unusually radio luminous stars; the sgB[e] star 
W9 and the RSG W26. In both cases the emission was found 
to be spatially resolved, suggestive of recent mass-loss events. 
Motivated by these unusual radio properties and the possibility 
of detecting emission from massive post-MS stars over a broad 
range of evolutionary stages we carried out a more extensive ra- 
dio observation of Wd 1 at four frequencies. 

We present the results of this survey in this paper. The obser- 
vations are described in section 2, with the radio sources iden- 
tified in section 3. Section 4 discusses the nature of the radio 
sources, section 5 is a brief discussion of the extended emission, 
and the results are summarized in section 6. 

2. Observations 

Observation of Wd 1 have been obtained at 8640, 4800, 2496 

and 1 384 MHz (corresponding to 3, 6, 13 and 20cm respectively) 
using the Australia Telescope Compact Array. A first epoch was 
obtained 1998 March 4-5 in the 6B configuration, followed by 
observations 2001 January 7-8 in the more compact 750C con- 
figuration and on 2002 May 18 (only 8640 and 4800 MHz) in 
the 6A configuration. Together, these observations ensure good 
coverage across the spatial frequency range observed, though do 
not have identical spatial frequency coverage at each frequency. 
Hence, caution is suggested when comparing flux measures at 
the different frequencies for resolved, extended emission. 

Observations of Wd 1 were interleaved with repeated obser- 
vations of the nearby bright point source 1657-56 in order to es- 
tablish the antenna-based complex gains. Initial editing and cali- 
bration of all observations was done using the Miriad data reduc- 
tion package (Sault et al. 1995). The gain solutions established 
using 1657-56 were subsequently "referenced" by interpolation 
to the Wd 1 observations. The absolute flux scale was estabhshed 
using the primary calibrator 1934-638, assuming fluxes of 2.84, 
5.83, 11.14 and 14.94 Jy at 8.6, 4.8, 2.5 and 1.4 GHz respec- 
tively. The fluxes derived for 1657-56 are shown in Table 1. No 
significant flux variations (> 10% of the source flux) were noted 
in preUminary analysis of the fluxes of the brightest unresolved 
sources (> 2 mJy at 8.6 GHz) in data from each epoch. Hence, 
data from the three epochs were combined into one dataset in or- 
der to improve the signal-to-noise ratio of the data and improve 
our ability to detect weaker sources. 

After phase-referencing and subsequent combination of the 
data from each epoch, deconvolution was done through visibil- 
ity model fitting, a technique widely used in VLBI image con- 
struction, rather than the commonly used clean technique. This 
was done using an automated modelling fitting routine, the mod- 
coNS macro within the smerf patch (Reid 2006) to the DmviAP 
package (Shepherd 1997). The major advantage of the "smear 
fitting" method implemented by smerf is that it yields higher 
resolution for significantly detected features than estimated by 



Table 1. Fluxes determined for the phase-reference 1657-56 





Config. 


8.6 GHz 


4.8 GHz 


2.5 GHz 


1.4 GHz 






(Jy) 


(Jy) 


(Jy) 


(Jy) 


1998 Mar 


6B 


1.54 


1.81 


1.92 


1.89 


2001 Jun 


750C 


1.44 


1.61 


1.71 


1.88 


2002 May 


6A 


1.63 


1.79 







the oft-used Rayleigh criteria (Reid 2006), represented by the 
FWHM of the synthesized beam. A model was established for 
the 8.6-GHz visibilities consisting of point and Gaussian source 
models, using the image data to guide the modelling process. 
This gave a model at the highest observed frequency, and hence 
at the highest resolution. Once the positions of the model compo- 
nents had been determined, the relatively high signal-to-noise ra- 
tio of the brightest sources in the field permitted the use of phase- 
only self-calibration to improve the antenna gain solutions from 
those derived from phase referencing alone. A "best" model was 
established with a final round of model fitting, keeping the posi- 
tions fixed, to refine the flux of the model components. Note that 
no amplitude self-calibration was carried out using these mod- 
els. This "best" model at 8.6 GHz was then used as the input 
model at the other frequencies. The resulting images at the four 
passbands are shown in Fig. 1 . 

3. The radio sources in Wd 1 

3.1. Identification of radio stars 

In order to identify stellar counterparts to the radio sources in 
Wd 1, an /?-band image (655 nm effective wavelength) of the 
cluster was obtained on 2004 June 1 3 using the FORS2 camera 
on the VLT. After re-projecting the optical image to the same 
projection as the ATCA images using the aips task ohgeo, the 
absolute position of the optical image was fixed by assuming 
that the position of the point source component of the 8.6-GHz 
emission from W9, the brightest source in the radio image, is co- 
incident with the peak of its optical emission, as determined us- 
ing AiPS routine maxftt. Since the position of 8.6-GHz emission 
components was established by phase-referencing, the absolute 
position accuracy of the model components is determined by the 
accuracy of the position of the phase-reference source 1657-56, 
and on the residual antenna gain phase as a result of the offset 
between the position of the phase-reference source and Wd 1, 
separated by 8.8 degrees. It is estimated this leads to an uncer- 
tainty of ~ 100 mas in the absolute positions of the radio sources 
(see Table 2). 

The resulting overlay of the 8.6-GHz image and the R-hand 
image is shown in Fig 2. A large number of optical sources are 
coincident with radio sources. Setting a 4(t point source detec- 
tion limit of 0.23, 0.43, 0.86 and 1.45 mJy beam ^ at 8.6, 4.8, 2.5 
and 1.4 GHz respectively, we identify the optical sources with 
associated radio emission in Table 2. 

The position of putative optical counterparts to the radio 
sources were also determined using the aips routine maxfet, and 
an offset between the radio and optical positions derived. The 
uncertainty in the relative position of the radio and the optical 
image is governed by how well the optical image is referenced 
to the radio image. By inspection of the position of W9, we 
estimate this to be an uncertainty of approximately ±200 mas. 
This is the dominant position uncertainty for the sources. Radio 
sources that have a position offset from a putative optical "coun- 
terpart" greater than 600 mas may not be associated, and so these 
potential mis-identifications are identified in Table 2. 
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Fig.l. The ATCA observations of Wd 1 at 8.6 GHz (upper left), 4.8 GHz (upper right), 2.5 GHz (lower left) and 1.4 GHz (lower right). In 
each image, contour levels are -3,3,6, 12,24,48,96, 192cr, where the Icr rms uncertainty is 0.06,0.1,0.2,0.4 mjy beam"' respectively at the four 
frequencies. These images were reconstructed using the nominal FWHM of the synthesized beam at each frequency, shown in the lower left-hand 
comer of each image. Note that these images are not corrected for the primary beam pattern. 



A total of 18 stars are identified as radio emitters, with 
optical-radio positional offsets < 0.6". We identify two ob- 
jects with associated or potentially associated radio emission for 
which no designation is available in the literature. Spectroscopy 
that could identify if they are bona-fide cluster members is not 
available for these objects so we do not assign Westerlund num- 
bers to them, but name them as D09-R1 and D09-R2 following 
the convention used by (Clark et al. 2008, hereafter C08). Three 
stars are identified as possible radio emitters (W5, W16a and 



D09-R2) since they are offset from the potential radio "coun- 
terpart" by between 0.6 and 1 .0 arcseconds, and some caution 
may be warranted identifying these stars with the nearby radio 
emission. 

With this many stellar radio sources, Wd 1 has the richest 
population of radio emitting stars known for any young mas- 
sive galactic cluster (e.g. Lang 2003; Moffat et al. 2002; Setia 
Gunawan et al. 2003). The stellar radio sources are blue, yellow 
or red super- or hypergiants and WR stars, representative of dif- 



4 



S.M. Dougherty et al.: Radio emission from the massive stars in the Galactic Super Star Cluster Westerlund 1 




Right Ascension (J2000) 

Fig. 2. 8.6-GHz image overlaid on a FORS i?-band image. The limiting magnitude of the /?-band image is » 17.5 magnitudes. The radio data are 
represented by contours with levels at -3,3,6, 12,24,48,96, 192cr, where the Icr rms uncertainty is 0.06 mjy beam"' . The radio sources with putative 
optical counterparts listed in Table 2 are identified by the circles and Westerlund numbers. 



ferent stages of massive star evolution. The supergiant B[e] star 
W9 is by far the brightest stellar radio emitter in the cluster, as 
anticipated from the "snapshot" observation of the cluster pre- 
sented by C98. Among the other four known blue hypergiants, 
only the LBV W243 is detected. Wd 1 has one of the richest 
populations of YHGs known in the Milky Way, with six mem- 
bers of this group (W4, 8a, 12a, 16a, 32, and 265), all of which 
are detected except for W8a and possibly W16a. All four known 
RSGs, six of the 16 WR stars in the field-of-view (24 known 
currently in the cluster), and four OB supergiants, are also asso- 
ciated with radio emission. 

A number of sources appear to consist of multiple compo- 
nents (W4, W9, D09-R1 and W237) with a central unresolved 
source surrounded by a more extended component that is of- 



ten significantly larger than the synthesized beam of the array. 
For W4 and W9, the extended component is clearly centred on 
the stellar source and an association with the extended emission 
is strongly implied. For D09-R1, the relationship of the stellar 
source and the extended emission is not as clear since the ex- 
tended emission is off'set from the associated point source. 

3.2. Radio star characteristics 

The radio fluxes of the sources were deduced from the model 
components derived in the model fitting process, which were 
then corrected for attenuation by the primary beam of the 
ATCA, dependent on their distance from the pointing centre. 
The FWHM of the ATCA primary beam is 5.86, 10.06, 19.91 
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Table 2. Characteristics of stellar radio sources in the Wd 1 cluster 



Source ID 


Spectral 


8.6 GHz'' 


4.8 GHz 


2.5 GHz 


1 .4 GHz 


Spec. Index 


Radio pos.*" 




Notes 


a(b) 


Type^ 


(mJy) 


(mJy) 


(mJy) 


(mJy) 




16 47 


-45° 


(") 




4(H) 
4 


F2 la+ 


0.8 ± 0.08 
2.2 ± 0.2r'' 


0.6 ±0.1 
2.8 ± 0.3r 


3.3 ± 0.3« 


3.8 ± 0.4 


+0.49 ± 0.31 
-0.23 ± 0.07 


1.445 


50 37.66 


0.37 


1 


9(C) 


sgB[e] 


24.9 ± 2.5 


18.1 ± 1.8 


9.5 ± 0.9 


7.8 ± 0.8 


+0.68 ± 0.07 


4.136 


50 31.41 






9 


30.5 ± 3.0r 


28.4 ± 2.8r 


21.0 ± 2.1r 


24.5 ± 2.5r 


+0.16 ± 0.07 










12a 


A5 la+ 


2.9 ± 0.3r 


2.9 ±0.3 


3.5 ± 0.4 


3.0 ±0.4 


-0.06 ± 0.08 


2.151 


50 59.23 


0.61 


1 


15 


O9.5-B0.51a 


0.6 ± 0.06 


0.7 ±0.1 


< 0.6' 


< 1.1 


> 0.0 


6.648 


50 30.31 


0.54 




17 


O9.5-B0.51a 


1.7 ±0.2 


1.6 ±0.2 


3.2 ± 0.3 


2.9 ± 0.4 


-0.42 ±0.11 


6.183 


50 49.40 


0.50 




WRB 


WN7 


4.3 ± 0.4r 


4.8 ± 0.5r 


4.4 ± 0.4 


4.0 ± 0.4 


+0.04 ± 0.07 


5.345 


51 5.21 


0.21 


2 


20 


RSG 


3.8 ± 0.4r 


4.2 ± 0.4r 


4.1 ± 0.4 


4.8 ± 0.4 


-0.11 ±0.07 


4.679 


51 23.81 


0.33 


2 


26(A) 


RSG 


20.1 ±2.0r 


21.9 ± 2.2r 


18.4 ± 1.8r 


7 


+0.07 ±0.11 


5.370 


50 36.57 


0.33 


2 


31b(WR V) 


WN8 


0.4 ± 0.06 


0.4 ±0.1 


< 0.6 


< 1.1 


0.00 ± 0.50 


3.821 


50 39.22 


0.45 




32 


F5Ia+ 


0.4 ± 0.06 


0.4 ±0.1 


< 0.6 


< 1.1 


0.00 ± 0.50 


3.675 


50 43.70 


0.31 




44(WR L) 


WN9 


0.4 ± 0.06 


< 0.3 


< 0.6 


< 1.1 


> 0.5 


4.210 


51 7.34 


0.23 




D09-R1 
D09-R1 .' 


BSG 


0.7 ± 0.07 
6.5 ± 1.2r 


0.8 ±0.1 
6.2 ± 1.2r 


11.0± 1.2r 


17.2± 1.7r 


-0.23 ± 0.27 
-0.61 ±0.11 


9.075 


51 10.23 


0.31 


1,3 


72(WR A) 


WN7 


0.5 ± 0.06 


0.5 ± 0.1 


< 0.6 


< 1.1 


0.00 ± 0.40 


8.332 


50 45.58 


0.00 




75 


RSG 


0.3 ± 0.06 


< 0.3 


< 0.6 


< 1.1 


> 0.0 


a.yii 


/in CO no 


0.46 




Zi /(rS) 


RSG 


1 o _i_ n o 
1.0 ± U.z 


i.o ± U.Z 


1 o _L n T 


T n _i_ n /I 
z.U ± U.4 


— U.Uu ± U.IU 


3.117 


52 18.99 


U.Zo 




237 


5.6 ± 2.2r 


7.5±2.1r 


10.6 ± 3.0r 


11.2± 1.9r 


-0.35 ±0.19 








1,3 


239(WR F) 


WC9d 


0.3 ± 0.06 


<0.3 


<0.6 


< 1.1 


>0.0 


5.196 


52 25.63 


0.50 




243(G) 


LBV 


1.5 ±0.2 


0.9 ±0.1 


<0.6 


< 1.1 


+0.87 ± 0.30 


7.496 


52 29.28 


0.32 




265(1) 


F5 la+ 


2.3 ± 0.3r 


2.1 ±0.2 


1.7 ±0.2 


1.8 ±0.4 


+0.20 ±0.11 


6.268 


49 24.19 


0.30 




5(WR S) 


WNlO-11 


0.3 ± 0.06 


< 0.3 


<0.6 


< 1.1 


> 0.0 


2.877 


50 19.76 


0.91 


4 


16a(L) 


A2 la+ 


1.6 ±0.3 


2.0 ± 0.2 


2.8 ± 0.3 


3.0 ±0.4 


-0.33 ± 0.09 


6.688 


50 41.99 


0.93 


4 


D09-R2 


BSG 


0.7 ± 0.06 


0.8 ±0.1 


1.2 ±0.2 


1.6 ±0.4 


-0.43 ±0.11 


6.867 


50 37.17 


0.70 


4 



a) Source numbers from Westerlund (1987). 

b) Alternate names from C02, Borgman et al. (1970),Negueruela & Clark (2005), and Crowther et al. (2006) 

c) Spectral types from C02, C05, Negueruela & Clark (2005) and Crowther et al. (2006) 

d) Fluxes derived from visibility modelling, corrected for the primary beam of the ATCA. 

e) J2000 positions of model components at 8.6 GHz, assuming reference source 1657-56 at a = 17:01:44.853, 6 = -56:21:55.964. The uncertainty 

in absolute radio position is estimated to be ±0.1 arcsec 

f) Separation of the optical and radio positions. 

g) For sources with two components at higher frequencies, when there is only one model component determined this represents a total flux at this 

frequency 

h) The symbol "r" denotes components that appear resolved by the ATCA at that frequency. 

i) Quoted upper limits are 3cr. 

1) Extended component 

2) These sources are embedded in extended emission which may influence the flux estimates. 

3) Not clear if the extended emission is related to the compact radio source or the underlying optical source. 

4) The radio emission is offset from the stellar position by more than 3x the position uncertainty introduced by phase-transfer i.e. 0.6 arcsec. 
Hence, association of radio emission with the optical source is uncertain. 



and 34.61 arcminutes at 8.6, 4.8, 2.5 and 1.4 GHz respectively. 
Given the size of the field-of-view is ~ 4 arcminutes, this leads 
to flux corrections of up to 1 .4x in the outer regions of the field- 
of-view at 8.6 GHz. In an attempt to allow for systematic errors, 
the flux uncertainties quoted in Table 2 are taken from the max- 
imum of the following three values: the Icr rms uncertainty in 
the image, 10% of the source flux, or for resolved sources the 
rms uncertainty in the image multiplied by the solid angle of 
the source. This represents a conservative approach to estimat- 
ing the flux uncertainties, and hence to the estimated spectral in- 
dices. However, it presumes that the model components used to 
describe a source are not influenced by other emission from the 
region surrounding the source. This may be an additional com- 
plication in determining the fluxes of several sources (WR B, 
W20 and W26) since they are embedded in more extended re- 
gions of emission that may impact the accuracy of the derived 
fluxes. In addition to the source flux, it is noted which sources 
appear resolved from the visibility modelling directly, avoiding 
the impact of visibility weighting during imaging on the clean 



beam size and the subsequent hmitations of source modelhng in 
the image plane. 

It is noted that the flux values determined in this study are 
approximately 10% higher than the values previously reported 
by C98. It is not clear if the apparent increase in the flux of W9 
is due to intrinsic source variations or due to the absolute cali- 
bration scale, given that the uncertainty in the absolute flux cal- 
ibration is ~ 5 - 10%. We suspect the apparent variation is due 
to a combination of absolute calibration uncertainty and the fact 
the data used in this study extend to lower spatial frequencies 
than in the study of C98, and hence may recover flux that was 
"resolved out" in the earlier observations. We suggest a more 
accurate determination of total flux is presented here than in the 
previous work. 

The radio spectral index a, where 5y oc v", for each source 
was calculated by a weighted regression fit of a single power- 
law to the fluxes. Use of a single power-law presumes there is 
no curvature in the continuum spectra over the frequency range 
observed. Inspection of the spectra of those sources detected at 
all bands shows no compelling evidence of curvature, given the 
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uncertainties in the radiometry at each band. The resulting val- 
ues for spectral indices and uncertainties are given in Table 2 and 
displayed in Fig. 3. For sources with radiometry at all frequency 
bands, the resulting single power-law index values are less af- 
fected by potential systematic errors in the flux measures at each 
band, compared to indices derived across a smaller frequency 
range e.g. 8.6 to 4.8 GHz. For sources that are undetected at one 
or more bands, and upper limits quoted, only the detected fluxes 
were used to calculate spectral index, unless they helped further 
constrain the spectral index value e.g. W44. It should be noted 
that the fluxes for sources only detected at one or two bands are 
all low (< 1 mJy) and the uncertainty in the derived spectral 
indices is high. 

None of the sources have a completely optically-thick ther- 
mal spectrum (a ~ 2), though several have spectra consistent 
with a ~ +0.6, expected for a partiaUy optically-thick, steady- 
state stellar wind e.g. the compact components in W4, 9, 44 and 
243. 

The majority of the sources have indices that are quite shal- 
low compared to the expected value for a stellar wind, with many 
having indices consistent with -0.1, the value for optically- 
thin thermal emission. A combination of optically thick and 
thin emission components can give rise to such spectra, such 
as clumpy stellar winds where the clumps are opticaUy thick 
(e.g. Ignace & Churchwell 2004). Alternatively, in some of these 
objects the emission may be due to a combination of thermal 
and non-thermal emission forming a composite spectrum, as 
observed in a number of massive stars, where the non-thermal 
emission is attributed to a coUiding-wind binary (CWB) (e.g. 
Chapman et al. 1999; Dougherty & WiUiams 2000; van Loo 
et al. 2006). Depending on the relative strength of the two contin- 
uum components, the spectral index could lie anywhere between 
— hO.8 and -0.6, this latter value typical of optically-thin syn- 
chrotron emission. This makes an unequivocal identification of 
the underlying emission difficult. To compound matters, com- 
posite spectra due to thermal-t-non-thermal emission may not 
be well represented by a single power-law (e.g. Pittard et al. 
2006). However, identification of CWBs is helped through ob- 
servations at other wavelengths, largely thermal X-ray emission 
arising from shock-heated plasma in the wind-collision region, 
and/or dust emission that appears to be a common feature in 
carbon-type WR star massive binaries with non-thermal emis- 
sion (Williams 1999). 

There are a few sources with indices that, at face value, are 
more negative than -0.1, suggesting non-thermal emission e.g. 
W17 and the extended component associated with W4. Other 
potential members of this group are W16a, D09-R2 and the ex- 
tended emission near D09-R1 but, as noted earUer, the associa- 
tion of the radio emission with the underlying stars is unclear. 
The derived spectral index for the extended emission in W4 is 
consistent with optically- thin thermal emission at the 95% con- 
fidence level. This leaves W17 as the only bona-fida non-thermal 
source. 

Additional observations over a broader range of frequencies 
are required to establish more firmly the nature of the radio emis- 
sion from the stars in Wd 1. Individual cases are discussed in 
Sec. 4 where the observations from each stellar sub-type are dis- 
cussed in detail. 

3.3. Extended radio emission 

In addition to the radio emission that is associated with stars it 
is clear there are a number of large, extended emission regions, 
distributed within ~ 1.5 arcminutes of the cluster centre. Some 



Table 3. Number of radio emitters of given spectral type 



Spec. 


No. radio 


Source 


Cluster 


Type 


emitters 


ID 


total 


OB SGs 


3(4?)° 


15, 17, D09-R1, (D09-R2?) 


-150 


sgB[e] 


1 


9 


1 


BHGs 


1 


243 


4 


YHGs 


4(5 ?r 


4, 12a, 32, 265 (16a?) 


6 


RSGs 


4 


20, 26, 75, 237 


4 


WN9-10 


1 (2?)" 


44, (5) 


2 


WN5-8 


3 


WRB. 31b. 72 


14 


WC 


1 


239 


8 



Notes: a) These sources are included in parentheses due to the uncertain 
association of radio emission with W16a, D09-R2 and W5. b) WRs N 
(WC8), O and Q (both WNL) are not in the fleld-of-view of the radio 
observations described here. 
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Fig. 3. The spectral indices of the stellar radio sources in Wd 1, iden- 
tified by their Westerlund numbers. The canonical spectral indices for 
partially optically thick thermal emission from a stellar wind (+0.6) and 
optically thin thermal emission (-0. 1) are highhghted by the dashed and 
dotted Unes. Error bars are 1 sigma. The arrows denote lower Umits to 
the derived spectral index. Radio sources offset from potential optical 
counterparts by more than 0.6" are marked by open boxes. 



Table 4. Properties of the extended radio emission regions in the cluster 
centre. 



ID 


8.6 GHz 


4.8 GHz 


2.5 GHz 


Spec. Index 


Assoc. 




(mJy) 


(mJy) 


(mJy) 




Star 


Al 


5.6 


7.6 


7.6 


-0.24 ±0.15 


D09-R1 


A2 


23.8 


26.0 


35.5 


-0.32 ±0.10 


WRB 


A3 


12.1 


11.1 


16.0 


-0.23 ± 0.20 


W20 


A4 


3.5 


4.7 


5.5 


-0.36 ±0.10 




A5 


2.0 


1.9 


2.0 


-0.01 ±0.11 




A6 


71.9 


78.0 


81.8 


-0.05 ± 0.07 




A7 


22.4 


28.2 


27.0 


-0.14 + 0.13 


W9 


A8 


153.3 


164.9 


161.4 


-0.04 + 0.11 


W26 


Total 


295 


322 


337 







Note: These fluxes are for the extended emission only in these regions. 
The fluxes of the radio stars in these regions has been subtracted. The 
flux uncertainty is assumed to be - 10%. 



of these regions appear to be associated with stars e.g. A3 with 
W20 and A8 with W26, but there are others such as A5 and A6 
for which optical counterparts are not identified in the FORS2 
image. These regions are identified in Fig. 4 and the radio prop- 
erties summarised in Table 4. The spatial frequency coverage is 
not the same at all frequencies and some extended emission may 
be resolved out at 8.6 GHz and hence the reported fluxes may be 
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low relative to the fluxes at lower frequencies. Hence, some of 
the spectral indices reported for the individual regions in Table 4 
may actually be flatter, or less negative, than reported. 

The total radio flux from Wd 1 as measured in the inter- 
ferometer data presented here is 422, 461, 523 and 669 mJy 
at 8.6, 4.8, 2.2, and 1.4 GHz respectively. This is compared to 
an interpolation of the single-dish fluxes reported by Kothes & 
Dougherty (2007), which gives fluxes at the same frequencies of 
450, 499, 560, and 620 mJy. Given that the flux uncertainties in 
the single dish measures is likely ~10%, the total fluxes derived 
here from the ATCA data are in excellent agreement with the 
single dish fluxes. 

To determine the flux in extended emission, we subtract the 
radio emission associated with the stellar sources as given in 
Table 2 from the total flux determined by the ATCA. This leads 
to fluxes of 307, 351 and 426 mJy at 8.6, 4.8 and 2.2 GHz, 
with an uncertainty of ~ 5 - 10% . The flux at 1.4 GHz is not 
quoted since it is difficult to ascertain accurately for the radio 
stars, most especially for sources embedded in extended emis- 
sion. Comparing these values to those given in Table 4 and tak- 
ing into account their uncertainty, shows good agreement, per- 
haps degrading toward 2.2 GHz, which suggests there is little 
other extended emission in Wd 1 beyond the regions demarked 
in Fig. 4. 

Using the total fluxes above for the extended emission 
in Wd 1 gives a spectral index of -0.26 + 0.07, consistent 
with optically-thin thermal emission. Following the prescription 
given in Sec. 4.3, an 8.6-GHz flux of 310 mJy from a region 1.5' 
in diameter leads to an ionized mass of 15 for a plasma tem- 
perature of 10 kK, in excellent agreement with the ionized mass 
derived by Kothes & Dougherty (2007). 

C98 showed that the nebular radio emission associated with 
the RSG W26 correlated with mid-lR emission. Figure 4 shows 
an overlay of the radio emission in Wd 1 with Spitzer GLIMPSE 
data at 8;um (Benjamin et al. 2003). Without being distracted by 
the saturation and artefacts in the mid-lR image, it is striking 
that all the extended radio emission is associated with mid-lR 
emission. Additionally, note that mid-IR emission is associated 
with each of the RSG and YHG stars in the cluster (Sec 4.3). 
Such a close correlation between the ionized emission evidenced 
in the radio and the 8//m emission is well established in ultra- 
compact Hii regions, where the 8^(m emission is due to excited 
PAH emission (e.g. Hoare et al. 2007) whereas in more evolved 
Hii regions there is a clear association with hot dust, and the PAH 
emission associated with the photo-dissociation region beyond 
the ionized gas (e.g. Povich et al. 2007; Watson et al. 2008) 

4. The stellar sources 

The post-MS evolution of stars with masses in excess of ~3OM0 
is characterized by significant mass loss. While such mass loss 
plays a critical role in determining the subsequent evolution 
of the stars, the physical processes driving it are currently ill- 
constrained, with the mass-loss rate and wind velocities antic- 
ipated to differ by several orders of magnitude depending on 
the exact evolutionary phase. Given that Wd 1 is uniquely well 
stocked with examples of every known post-MS stellar type, 
these observations offer the possibility of constraining this pro- 
cess. 

Consequently, the following discussion on the nature of 
the radio emission for the 18 stellar detections is presented in 
a likely evolutionary sequence from OB supergiants, through 
cool super-/hypergiants, early hypergiants and WN and WC 
Wolf-Rayet stars. Given their comparatively weak winds (~ 
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Fig. 4. Comparison of radio emission at 8.6 (green) and 4.8-GHz (red) 
observations with 8/jm emission taken as part of the Spitzer GLIMPSE 
project. Only a few radio contours that demark the extended radio emit- 
ting regions are shown for clarity. 

10"^ Mq yr"'), the mid-to-late O-type main sequence progen- 
itors of these stars are not detected (Martins et al. 2005). The 
emission associated with the sgB[e] star W9 is discussed out- 
side this evolutionary sequence, due to the uncertain placement 
of the sgB[e] phase in the post-MS stellar zoo. 

Throughout the discussion it will be assumed that Wd 1 is at 
a distance of 4 kpc. Recent deep-lR imaging identifies the MS 
and pre-MS populations of Wd 1, leading to a photometric dis- 
tance of 4.0+0.3 kpc (Brandneret al. 2008), consistent with both 
4.7 + 1.1 kpc from an analysis of near-lR photometry of the WR 
stars in Wd 1 (Crowther et al. 2006), and on an atomic hydrogen 
absorption distance of 3.9 ±0.7 kpc (Kothes & Dougherty 2007). 

4.1. W9 - a luminous radio source 

By far the brightest radio source in Wd 1 is W9, with a to- 
tal flux at 8.6 GHz of 55.4 mJy, giving an 8.6-GHz luminos- 
ity of 1.6 X 10^' erg s"'. This places W9 among the most ra- 
dio luminous thermal emitting stars known: a factor of a few 
less than the 8.6-GHz luminosity of the extreme LBV 77 Carina 
(4.5 X 10^' erg s"') at radio minimum (~ 600 mJy; Duncan & 
White 2002), comparable to the massive overcontact binary RY 
Scuti (Gehrz et al. 1995) and a factor of a few more luminous 
than the other oft-cited radio luminous thermal emitter MWC 
349A (e.g. Tafoya et al. 2004). 

C98 hypothesised that W9 comprised of two components, a 
compact stellar wind source and an extended region. The obser- 
vations presented in this paper support this model, with an un- 
resolved source surrounded by an extended region of emission 
that are both detected at all four observing frequencies. 

The spectral index of the unresolved source coincident with 
W9 is +0.68 + 0.08, consistent with thermal emission arising 
from a partially optically thick, steady-state stellar wind with 
an r^^ radial density distribution. It is assumed this component 
represents the present day stellar wind from W9. The spectral 
index of the extended region is +0.16 ± 0.07, which we take 
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Table 5. Best-fit parameters of the r ^ shell model applied to the con- 
tinuum spectra of the compact and extended radio components of W9. 



Model 


Component 


M Meyr"' 
(xlO-5) 


(AU) 


K 

(AU) 




1 
2 


Compact 
Extended 


^•'-0.2 

17 6+"-^ 

^ '-"-0.4 




iioo!^« 


5800!-S 

00 


1 


3 
3 


Compact 
Extended 


Q 9+U.4 

33:^^ 




4QOO+3000 
"""-■win 


^"""-900 
OO 


3 
IL 



Notes: Models 1 and 2 correspond to best-fits of each individual spectra. 
Model 3 assumes the outer radius of the compact component is equal 
to the inner radius of the extended component. For each model it is 
assumed that T, = 10 kK, v„ = 200 kms"',Z = 0.9,ju = 1.5, y = 0.8 
and the volume filling factor / = 1 . 



to be essentially flat, and arguably consistent with optically-thin 
thermal emission. Assuming the extended region has a r^^ radial 
ion density distribution, similar to the stellar wind, the lack of a 
turnover in its continuum spectrum to a positive spectral index 
impUes that the extended region is totally optically thick down to 
1.4 GHz. With the free-free opacity for an ion distribution at 
frequency v behaving as t,,(/7) oc /j^-'v"^ ' where p is the impact 
parameter (see Appendix A and Eqn. A. 2), the lack of a turnover 
impUes that the extended region must have an inner radius that 
is larger than the radius of the Ty = 1 surface at 1.4 GHz. In this 
case, it is suggested that the extended region represents an earlier 
phase of mass-loss from W9, prior to the start of the current 
stellar wind phase. 

To model the continuum spectra of both the stellar wind and 
the extended envelope a shell-like geometry is adopted, with 
outer and inner radii Rg and respectively and with a r^^ radial 
ion distribution (e.g. Taylor et al. 1987). The free-free opacity at 
radius p in such an envelope is determined using Eqns. A.2, A.3 
and A.4. Together with Eqn A.l, these lead to the flux from 
the envelope at frequency v. Assuming an electron temperature 
Te = 10 kK and an outflow velocity of 200 kms a best-fit of 
this model to the continuum spectra of the two components over 
the four observing frequencies gives the parameters of the two 
components (Table 5 and Fig. 5). The uncertainties in the fitting 
parameters are derived in the standard manner of fixing all but 
one parameter which is varied until changes by unity. 

A total of three models were examined. Model 1 leads to 
a mass-loss rate of 8.7+°-^ x 10-5(Voo/200 km&-^)f/^ Moyr"! 
forZ = 0.9,a; = 1.5, r = 0.8,r = 10 kK (e.g. Leiflierer & 
Robert 1991) for a volume filling factor for clumps / (Abbott 
et al. 1981). The mass-loss rate is scaled to 200 km s ' since the 
outflow velocity in W9 is not well constrained with observed Une 
widths that range between 40 and 800 km s"^ (COS). Assuming 
the same outflow velocity and clump structure, the derived mass- 
loss rate of the extended region (Model 2) is approximately twice 
as high as the current phase of mass-loss, implying an earher era 
of higher mass-loss rate ~ 10"^ Moyr"^ 

The model spectrum of the individual radio components es- 
tablished that the best-fit outer radius of the compact compo- 
nent (the current stellar wind) is greater than the inner bound 
of the extended region (earher phase of mass-loss). If the mass- 
loss that produced both components is radially symmetric, it is 
not plausible for the inner radius of the extended region to be 
smaller than the extent of the current stellar wind. The third 
model (Model 3) assumed the inner bound of the extended re- 
gion is at least as large as the outer bound of the current stel- 
lar wind. Fitting both continuum spectra simultaneously through 
minimizing the total from the two models gives the parame- 
ters of Model 3, which supports a factor of 3 higher mass-loss 



Frequency (GHz) 



Frequency (GHz) 

Fig. 5. Model fits to the continuum spectra of the unresolved partially- 
thick stellar wind (top) and extended, optically-thin (bottom) compo- 
nents of W9. The solid lines represent the best-fit model when the stel- 
lar wind outer radius/extended inner radii are constrained to be identical 
(Model 3). The dotted spectra are the best-fit models when the stellar 
wind outer radius and extended region inner radius are free parameters 
(Models 1 and 2). 



rate in the earlier phase of mass-loss than in the current stellar 
wind. The resulting model fits to the data are shown in Fig. 5. 

The derived mass-loss rate for the putative earher mass-loss 
phase is high for a massive star, being of the same order as the 
the current mass-loss rate in t] Carina of » 10""* - 10"^ M© yr"' 
(Corcoran 2008). The inner bound of ~ 4000 AU suggests that 
this mass-loss epoch would have ended ~ 200 yrs ago assuming 
an outflow velocity of 200 km s"* , followed subsequently by the 
current mass-loss phase with 3x lower density. 

With such a model for the circumstellar environment, what 
does this imply for the evolutionary status of W9? While veil- 
ing any photospheric features, the rich optical-lR emission line 
spectrum and IR excess indicate a rich circumstellar environ- 
ment and permit a formal sgB[e] star classification (C08; Ritchie 
et al. 20 1 0). However spectral peculiarities exist when compared 
to other examples; most notably the unusual composite line pro- 
files of He 1 5876 and 6678A. The X-ray emission from W9 sup- 
ports a binary system, being too hard (kT ~ 3 keV) and bright 
(Lf ~ 10^^ ergs~^) to arise in a single star, although it does 
not constrain the nature of the components (C08; Ritchie et al. 
2010). Likewise, while W9 is found to be photometrically vari- 
able, no period has been identified currently (Bonanos 2007). 
Identification of W9 as a binary awaits further supporting obser- 
vational evidence. 



S.M. Dougherty et al.: Radio emission from the massive stars in the Galactic Super Star Cluster Westerlund 1 



9 



The current mass-loss rate of W9 is approximately an order 
of magnitude higher than that inferred for any of the other tran- 
sitional stars within Wd 1. It is directly comparable to those of 
the LBVs AFGL 2298, AG Car, FMM 362 and the Pistol star 
in quiescence with M = 3 - 6 x 10"^/'^^ Moyr"' (Clark et al. 
2009; Groh et al. 2009; Najarro et al. 2009), despite being a fac- 
tor of ~ 2x less luminous'. Moreover, depending on the outflow 
velocity and wind-clumping factor the mass-loss rate for W9 is 
uncomfortably close to the limit expected for a line-driven wind 
for which Smith & Owocki (2006) estimate a maximum mass- 
loss rate of 

1.4 X lO-\L/lO^Lo) Moyr"'. 

Such a high mass-loss rate and the presence of a detached 
shell of material formed in a previous phase of enhanced mass 
loss is suggestive of an LBV identification, with M ~ 3 x 
10"^ Mgyr"' for W9 during this event matching the time aver- 
aged mass-loss rate inferred for the LBV Wra751 by Voors et al. 
(2000), as weU as the YHG IRC-i-10 420 during the formation of 
its nebula (Sect. 4.3). The observed line spectrum, absolute vi- 
sual magnitude and spectral energy distribution exclude the pos- 
sibility that W9 is currently in a cool hypergiant phase with the 
binary companion providing the requisite ionizing flux to yield 
the emission line spectrum. Therefore, if the enhanced mass-loss 
rate is solely driven by the underlying star, this would be one of 
the few known examples of an LBV caught in a high mass-loss 
rate episode. Moreover, it would demonstrate that significantly 
enhanced mass loss in ~ 4OM0 stars is not solely Umited to oc- 
curring in a cool hypergiant phase. 

However, while the current observational data are consistent 
with the ejecta being formed as a result of the post-MS evolution 
of a single massive star, its formation in an interacting binary 
system undergoing significant mass transfer, in turn leading to 
common envelope evolution, also appears viable and well mo- 
tivated given the evidence of W9 being a binary. Examples of 
massive binaries with both a normal and a relativistic compan- 
ion are known: RY Scuti and SS433 respectively, and thus pro- 
vide observational templates. As with W9, both these systems 
have rich emission-line spectra and IR excesses due to the pres- 
ence of circumstellar dust (Gehrz et al. 1995; Clark et al. 2007, 
and references therein). More interestingly, both are associated 
with compact, bright radio nebulae with sizes and fluxes directly 
comparable to W9 yet attributed to mass lost through binary in- 
teraction (BlundeU et al. 2001; Gehrz et al. 1995). The former 
authors derive a mass-loss rate of the order of 10""' Mgyr"' from 
the radio observations for SS433, while the latter infer an ionized 
mass of lO"-' Mq for the nebula associated with RY Scuti, with 
a mass-loss rate of ~10"'* - lO""" Mgyr"' inferred from optical 
emission fines (de Martino et al. 1992); both rates are compara- 
ble with W9. 

4.2. The OB supergiants 

Currently, there are -75 spectroscopically identified evolved O 
stars (O giants and OB supergiants) within Wd 1, with a com- 
parable number of stars classified as such on the basis of pho- 
tometry (C05; COS; Negueruela et al. 2010)^. Unfortunately, the 

' In the absence of a bolometric correction for W9 adopting a compa- 
rable luminosity to the cool h3'pergiants - appropriate considering post- 
MS evolution is expected to proceed at constant bolometric evolution 
prior to the WR phase - suggests log(L/Lo)~5.8. 

^ Note that wc include implicitly the mid-B supergiants such as 
WTO, 71 and 57 here, while deferring discussion of the more lumi- 



currently published spectroscopy of these objects is of insuffi- 
cient resolution and wavelength coverage to yield accurate phys- 
ical parameters. However, canonical wind parameters for late 
O/early B supergiants (see C05, Table A2) indicate that at a dis- 
tance of ~ 4 kpc thermal wind emission would be undetectable, 
due to a combination of comparatively high terminal wind ve- 
locity and low mass-loss rate. As an example, the 09.5 la star 
with Voo ~ 1900 kms"' and M = 6 x 10"^ Mgyr"' and adopting 
Z - y = I and jj = 1.34 (Leitherer et al. 1995) has an expected 
8.6-GHz flux of ~ 0.06 mJy, much lower than the 8.6-GHz de- 
tection limit of 0.17 mJy. Since radio flux is oc f~^^^, it would 
take a highly clumpy wind with / ~ 0.5 to make such a wind 
detectable. 

Consequently, the apparent detections of W15 and W17 
(091b and 09Iab respectively - Negueruela et al., in prepara- 
tion), D09-R1 and potentially D09-R2 (both BOI or earlier; COS) 
are somewhat unexpected, most especially considering that these 
stars are unremarkable compared with the rest of the OB stars 
studies spectroscopically. Of note is that W15, W17 and D09- 
Rl all are coincident with extended emission regions and the pa- 
rameters of the unresolved radio components coincident with the 
underlying stars are likely influenced by this extended emission. 
As a result, the derived parameters may not reflect the intrinsic 
properties of the underlying stars, but rather the extended plasma 
region. 

That withstanding, W17 has a spectral index marginally 
steeper than -0.1, hinting at non-thermal emission and a poten- 
tial CWB. However, we note it does not have detected X-ray 
emission COS and shows no evidence of RV variations (Ritchie 
et al. 2009a). W15 has a spectral index lower limit of 0.0, and 
it is not possible to distinguish whether the emission is from ei- 
ther a stellar wind or a non-thermal/thermal composite source. 
It shows no evidence for binarity in either RV observations 
(Ritchie et al. 2009a) or X-ray observations (COS), where the 
emission is soft (kT ~ 0.6 keV) and at a level expected for a sin- 
gle O-type star (~ 2 x 10^^^ erg s '; e.g. Berghoefer et al. 1997) 
as a result of shock heating within the stellar wind. Nevertheless, 
if the radio flux is attributed to thermal emission from a stel- 
lar wind, W15 would have an unusually high mass-loss rate of 
M = 3 X 10"^ Moyr"' in comparison to other late-O/early-B 
supergiants (Crowther et al. 2006; Mokiem et al. 2007). The ab- 
sence of RV variations in these stars does not rule out a binary 
interpretation unequivocally since the mass ratios may be ex- 
treme and the inclination of the orbit unfavorable for ready de- 
tection of RV variations. Furthermore, the orbital period could 
be significantly in excess of the observation timescale. 

The spectral index of the point source associated with D09- 
Rl is comphcated by the presence of an extended component, 
and given the absence of both optical spectroscopy and an X-ray 
detection, it is simply concluded that the nature of this emission 
remains uncertain, and note that if it is assumed to be a stellar 
wind, the same problems with respect to an extreme mass-loss 
rate as found for W15 are encountered. 

4.3. The cool super- and hypergiants 

All four of flie RSGs and at least four of the six YHGs in Wd 1 
are radio sources. Irrespective of the uncertain identification of 
W16a as a radio emitter, the fraction of radio sources amongst 
both these classes of object is high. Of the RSGs, W75 was the 
weakest source detected and was unresolved at any wavelength. 



nous, and likely more evolved, B hypergiants W7, 13, 33 42a and 243 
to Section 4.4. 
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The other three detected RSGs (W20, 26 and 237) have extended 
radio emission. 

For both W20 and W26 the radio emission is best charac- 
terised by two components: a compact (~ 2"), spatially resolved 
source coincident with the star, and a more extended nebula 
(> 10"). W237 is similar in structure though the inner compo- 
nent coincident with the star is unresolved. In W20 and W26 
the more extended components have a distinct cometary mor- 
phology that extends away from the central region of the cluster, 
while that in W237 is more elliptical. 

The radio properties of the YHGs are more heterogeneous. 
W8 is not detected. W32 is associated with a weak point source, 
while W12 and W265 are associated with compact resolved 
sources, similar to the RSG radio sources. Finally, the radio 
emission from W4 consists of a point source at the position of 
the star surrounded by an extended component. 

The majority of the cool super- and hypergiants have spec- 
tral indices that are flat, within the uncertainties. Unlike OB 
and WR stars, composite thermal and non-thermal spectra have 
not been identified previously, and are not expected. Hence we 
identify these indices with optically-thin thermal emission. For 
a Maxwellian plasma at temperature T, the opticaUy-thin radio 
flux density is given by 

Sy = 5.7 X \Q-^^T^^^gyD-^Ev mJy, 

where D is the distance in kpc and Ey is the emission integral 

at 8.6 GHz. The total ionized mass can then be determined from 
Mt ^ luniHiEv/rie). 

The ionized masses of the nebulae, based on the source pa- 
rameters determined at 8.6 GHz, are given in Table 6. Note 
this is not the total mass since these observations are insensi- 
tive to any neutral component. On average, the emission integral 
Ey ~ 1.5 X 10^** cm"-', implies an ionizing luminosity that has 
to be greater than EycP'^ ~ 6 x 10^^ s ^ where a^^' is the is the 
hydrogen recombination rate to all but the ground state. Such 
ionizing luminosity is not available from these relatively cool 
stars with stellar efl'ective temperatures < 5000 K (e.g. de Jager 
1998). 

How are these envelopes ionized? A cluster of mass lO^M© 
and solar metalUcity can easily provide such a photon luminos- 
ity and seems the most likely source of the ionizing radiation. 
Taking the cluster to be approximately 1 pc across and with 
~ 100 OB-type stars implies that on average each OB star oc- 
cupies a volume 0.04 pc-', and is separated from the next OB star 
by ~ 0.34 pc. If each OB star provides 10"*^ s"' ionizing photons 
(Martins et al. 2005), with ~ lO'*^ s"^ required to fully ionize the 
cool hypergiant envelopes, suggests a radiation dilution factor 
oc (Rt/D)~^ less than 1000, which implies the hypergiant en- 
velopes are ~ 0.01 pc radius. This is closely consistent with the 
sizes observed for RSG nebulae e.g. VY CMa (see Sec. 4.3.2). 

Alternatively, a hot companion star could provide the re- 
quired photon flux, which at 10"*^"^^ s~^ implies an early B-type 
star. The YHG HR 8752 provides a precedent for such a model: 
it is detected at radio wavelengths and together with [Nn]- 
emission lines, needs a source of ionizing radiation that is pro- 
vided by a putative Bl companion (Stickland & Harmer 1978). 
In this respect we highlight the presence of [Nn] -emission lines 
in the optical spectra of W20, W237 and W265 (Clark et al. 
2010). 

If the plasma in these optically-thin nebulae is due to mass 
loss from the underlying stars, having optically-thin spectra 



Table 6. Ionized mass estimates for the optically-thin stellar nebulae 
around the YHGs (top panel) and RSGs (bottom panel) in Wd 1. 







Maj 


Mm 
Major 


VlJl 

10» 




10' 


iVlaSa 

10-' 




(mJy) 


(") 




(cm') 


(cm-') 


(cm-') 


(Mo) 


4 


2.2 


2.45 


0.49 


4.0 


1.3 


5.7 


2.9 


12a 


2.9 


1.82 


0.92 


5.7 


1.7 


5.5 


3.9 


265 


2.3 


2.86 


0.54 


7.7 


1.4 


4.2 


4.1 


20 


15.6 


13.79 


0.48 


677.9 


9.2 


1.2 


99.6 


26 


173.4 


15.0 


0.33 


412.4 


102.2 


5.0 


259.0 


237 


5.6 


11.19 


0.76 


908.0 


3.3 


0.6 


69.1 



Note: A distance of 4 kpc is assumed 



down to a frequency of 1.4 GHz imphes that the ionized ex- 
tended envelopes have inner radii significantly larger than the 
stellar radii. If this is the case, two possible models for the en- 
velope can be advanced. The mass loss from these objects could 
have stopped or significantly decreased at some point in the past 
resulting in a detached shell. Alternatively, the ionized region 
may only occupy the periphery of the envelope, with the inte- 
rior part of the envelope being neutral. Such a geometry could 
result from external ionization from the cluster radiation field. 
High-resolution observations are required to substantiate either 
of these potential models. 

4.3.1 . The Yellow Hypergiants 

Given the rarity of YHGs and the difficulty identifying suitable 

observational diagnostics, few measurements of the quiescent or 
outburst wind velocities and mass-loss rates of YHGs exist in 
the literature. For HD 33579 (A31a+) in the LMC, Stahl et al. 
(1991) quote Voc ~ 200 kms"' and M ~ 2 x 10"*^ Moyr-\ 
and Israelian et al. (1999) report Voo ~120kms"^ and itf < 
6.7 X lO-^Moyr ' for HR 8752 in 1998. For the well-known 
example p Gas, Lobel et al. (1998) find Voo ~ 120 kms-' and 
M < 9.2 X IQ-^Mq yr-' in 1993, in contrast to v„ ~ 35 km s-' 
and M ~ 5.4 X 10"^ Mgyr"' during the 2000 outburst (Lobel 
et al. 2003). Castro-Carrizo et al. (2007) measure outflow veloc- 
ities of 35 kms-' and 25 - 37 kms-' for the ejecta associated 
wifli HD179821 and IRC -hIO 420 (although a priori it is not 
clear if they were ejected in a RSG or YHG phase). For both 
these stars, variable mass-loss rates of 4 - 300 x 10"^ Mgyr"' 
and 1.2 - 3 x lO"'* Moyr ' were found respectively during the 
nebulae formation events, with Blocker et al. (1999) also report- 
ing a mass-loss rate in excess of 10-^ Moyr^' for IRC -hIO 420 
in the recent past. These works all suggest that during outburst 
the mass-loss rates in YHGs are higher than during quiescence 
whereas the outflow velocity is lower during outburst. 

Taking a mean 8.6-GHz flux of 2.5 mJy for the YHGs, their 
radio luminosity is 4.9 x 10'^ erg s"'. The only other radio de- 
tected YHG in the Galaxy, HR8752, has an 8.6-GHz flux of 
17 mJy (Dougherty et al., in preparation) at a distance of 1 .7 kpc, 
giving a luminosity of 5.9 x 10'^ erg s-', similar to the YHGs in 
Wd 1. HR 8752 has been interpreted to have a stellar wind from 
measured spectral indices of 0.51 ± 0.08 (Higgs et al. 1978) and 
0.61 + 0.04 (Dougherty et al., in preparation). Only the compact 
component in W4 appears to have a spectral index consistent 
with a stellar wind (h-0.49 + 0.31), but given the uncertainty, 
arguably also consistent with optically-thin thermal emission. 
Assuming Z = 0.9, ju =1.5 and y = 0.8 as for W9 and that 
the putative wind from W4 is fully ionized and smooth (/ = 1) 
gives a mass-loss rate of 0.7xl0-^(Voo/200km s"') Moyr-'. This 
is consistent within a factor of a few with the quiescent mass-loss 
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rates of HD 33579, HR 8752 and p Cas, but significantly lower 
than inferred for p Cas, HD179821 and IRC +10 420 during out- 
burst, especially for comparably low velocities (< 40km s~^). 

The extended nebulae around W4, W12a and W265 have 
masses of 3 - 4 x 10"^ M© and radii of ~ 0.02 pc (based on 
the major axis size in Table 6). An assumed outflow velocity of 
~ 200 km s"^ gives a flow time of ~ 100 yrs, during which time 
mass loss at a rate of ~ 10"^ Mgyr as found for W4 above, 
yields a total mass loss of 10"^ M©. Thus, by comparison with 
the cases of p Cas, IRC +10 420 and HR 8752, the radio nebulae 
associated with the YHGs in Wd 1 may be the result of quiescent 
mass loss, rather than an outburst phase. 

None of the YHGs in Wd 1 possess the old (~ 4000 yr), mas- 
sive (M> 1 Mq) ejection nebulae seen in HD179821 and IRC 
+ 10 420 (Oudmaijer et al. 1996; Castro-Carrizo et al. 2007). We 
suspect that in the extreme environment of Wd 1 such nebulae 
would quickly be ablated by the cluster wind arising from the 
OB stars - as appears to be the case for the RSGs W20 and W26 
(see Sec. 4.3.2). However we cannot exclude the possibility that 
the YHGs may be evolving to cooler temperatures (redwards) 
and have yet to encounter the instabilities leading to extreme 
mass-loss rates in either the RSG or post-RSG phase. A mixture 
of both pre- and post-RSG objects combined with the intrinsic 
variability of YHGs (Clark et al. 2010) may explain the hetero- 
geneous radio properties of the Wd 1 population. This argument 
is supported by evidence that even the best two candidates for 
post-RSG stars, HD 179821 and IRC +10 420, exhibit quite dis- 
similar nebular properties (Oudmaijer et al. 2009; Patel et al. 
2008). 

4.3.2. The Red Supergiants 

Given very high mass-loss rates inferred for extremely lumi- 
nous RSGs (log(L/Lo)> 5.5), it is suspected that they play an 
important role in post-MS evolution of stars at the Humphreys- 
Davidson limit. Consequently, they have been the subject of nu- 
merous (multi-wavelength) studies to constrain their mass-loss 
histories. Of the Galactic examples, Schuster et al. (2006) used 
tiie HST/WFPC2 to identify compact (r < 0.003pc) nebulae 
associated with NML Cyg, VX Sgr and S Per, with Morris & 
Jura (1983) also associating a much larger (>lpc) Hii region 
with NML Cyg. de Wit et al. (2008) report the presence a neb- 
ula around // Cep at 25 microns (r ~ 0.03 pc); a comparably 
sized nebulae surrounds VY CMa (r ~ 0.04 pc in HST/WFPC2 
images; Smith et al. 2001). Finally Serabyn et al. (1991) and 
Yusef-Zadeh & Morris (1991) report the presence of an asym- 
metric nebula (~lpc in length) around the Galactic Centre RSG 
IRS7 at both near-IR and radio wavelengths. 

As with these stars, the discovery of large nebulae associated 
with three of the 4 RSGs in Wd 1 (major axis > 0.2 pc; Table 6) 
emphasises the role of mass loss in this evolutionary stage. NML 
Cyg and IRS7 are of special interest in the context of the RSGs 
in Wd 1, given that these stars are also located within dense 
stellar aggregates, namely Cyg OB2 (Morris & Jura 1983) and 
the Galactic Centre cluster (e.g. Yusef-Zadeh & Morris 1991) 
respectively. Indeed the pronounced cometary nebulae of W20 
and W26 appear similar both in morphology and Unear scale to 
that of IRS7. Consequently we attribute the shape of these neb- 
ulae to a similar physical process, namely the interaction of the 
ejecta/RSG stellar wind with the cluster medium/wind. Lacking 
dynamical information, we refrain from attempting a quantita- 
tive analysis of the nebular morphologies, but note that the un- 
remarkable radial velocities found for these stars by Mengel & 
Tacconi-Garman (2007) suggest that the dominant cause of ram 



pressure shaping the nebulae is the expansion of the cluster wind, 
rather than the rapid motion of the stars through a quasi-static 
cluster medium. This is supported by the orientation of each of 
these cometary nebulae away from the cluster centre region. 

The morphology of the radio nebula associated with W237 
suggests it has been less affected by the cluster wind than either 
W20 or W26. Adopting an outflow velocity of 30 km s ' for the 
RSG wind, flie mass of 0.07 Mq and radius of ~ 0. 1 1 pc (Table 6) 
results in a kinematic age of 3,600 yr and a time averaged mass- 
loss rate of 2x10"^ Moyr"'. This age supports the suggestion 
that the radio nebula is little affected by the cluster wind, oth- 
erwise it Ukely would be ablated as observed in W20 and W26. 
This is consistent with its position further from the cluster centre 
than either W20 or W26. The derived mass-loss rate is compa- 
rable to the mass-loss rates found for field RSGs, which range 
from 10-* - 10"'' Moyr-i (Jura & Kleinmann 1990; Sylvester 
et al. 1998), although transient rates as high as ~10~^ Moyr"' 
(NML Cyg and VY Cma; (Blocker et al. 2001 ; Smith et al. 2001) 
and as low as a ~few 10"^ M© in Cep (de Wit et al. 2008) have 
been observed. 

Moreover, both W20 and W26 have similar nebular masses 
to W237 and, in turn, all three are directly comparable to 
that of the nebula surrounding VY CMa (Smith et al. 2001). 
Consequently the ejecta associated with the RSGs within Wd 1 
appear to be the result of similarly extreme, possibly transient 
mass-loss events as those that yield the nebulae around other 
Galactic RSGs (e.g. Smith et al. 2001). In this respect, the ap- 
parent instabilities in W26 and W237 are of interest (Clark et al. 
2010). 

4.4. The B hypergiants 

The final subset of transitional stars to be considered are the B 
hypergiants, which likely represent the immediate progenitors 
of WN-type stars. Included in this group are the LBV W243 and 
the cool B-type hypergiants W7, 33 and 42. The final member is 
the early-B hypergiant W13, which emphasises the close physi- 
cal continuity between the B hypergiants and the hydrogen-rich 
late-WN stars such as W44 (=WR L (WN9h); COS) detected at 
8.6GHz (Table 2 and Sec 4.5). Of these stars, only W243 has a 
spectral index indicative of thermal wind emission. 

W243 has recently been the subject of an extensive spectro- 
scopic study by Ritchie et al. (2009b). This confirmed the finding 
of C02 that W243 had substantially cooled since the previous ob- 
servations of (e.g. Westerlund 1987), with the stellar temperature 
between 2002-8 being 8-9kK, with a formal spectral classifica- 
tion of mid-A type and a photospheric absorption spectrum of 
a YHG. However, detailed non-LTE analysis were unable to si- 
multaneously reproduce these features and the prominent H i and 
He I emission lines, leading Ritchie et al. (2009b) to conclude 
that W243 harbored a hot binary companion. Unfortunately, cur- 
rent spectroscopic data are unable to constrain the nature of the 
hot companion; nevertheless at this time W243 closely resem- 
bles the known YHG binary HR8752 (Sec. 4.3.1) and we as- 
sume that the hot companion is responsible for the ionization of 
its wind. 

Given the current properties of W243 we adopt identical val- 
ues of Z, p and y as used for the YHG W4, leading to a mass-loss 
rate of 1.1 x 1 0"^ (v.^, /200km s"^)/'/^ Moyr"', which is directly 
comparable to W4. This is to be expected considering their simi- 
lar radio fluxes (Table 2). This mass-loss rate is an order of mag- 
nitude higher than that found by Ritchie et al. (2009b), correcting 
for difference in assumed outflow velocity, though it should be 
noted that the observed optical spectrum is relatively insensitive 
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to the current mass-loss rate. Assuming a 4£r-detection limit of 
0.24 mJy for the remaining stars yields a corresponding upper 
limit to the mass-loss rates of 2 x 10"^(Voo/200km s"') M^yr'^ 
In comparison, the LBVs HR Car and HD 160529 have 
mass-loss rates of ~ 1 - 2 x 10"^/'^^ Moyr ' determined from 
radio observations, with HD 80077 and HD 168607 an ~order of 
magnitude lower (Leitherer et al. 1995). Therefore, we conclude 
that the mass-loss behavior of the Wd 1 population of B hy- 
pergiants is consistent with expectations drawn from the wider 
Galactic population. No extended emission is associated with 
any of these stars, whereas a number of Galactic LBVs are as- 
sociated with radio bright nebulae. However such nebulae are 
typically found to be old (~ 10^ - lO'' yr Clark et al. 2003) and 
of an extent sufficient to encompass the central regions of Wd 1 . 
As argued for the YHGs, we suspect that they would quickly be 
disrupted by the cluster wind and/or radiation field. 

4.5. Wolf-Rayet stars 

In recent years, Wd 1 has been shown to harbour at least 24 
WR stars (Crowther et al. 2006), representing 8% of the cur- 
rently known galactic population. Five of these have been de- 
tected in this radio survey at a flux level typically 0.3 - 0.5 mJy 
at 8.6 GHz: WR A (=W72), B, F(=W239), L(=W44), and 
V(=W31b) (Crowther et al. 2006, and references therein), po- 
tentially consistent with the radio emission originating in their 
stellar winds. 

WR B is an exception, with a deduced flux of 4.3 mJy. This 
flux is 5 - lOx higher than the other detected WR stars, imply- 
ing a wind density that is 3 - 6x higher than the mean value for 
the other WRs. Bonanos (2007) found that WR B is an eclipsing 
binary and in binary systems the flux can be raised substantially 
above that from a stellar wind by emission arising in a wind- 
colUsion region (WCR) from the colliding winds of two mas- 
sive stars in a binary (e.g Dougherty et al. 2005). More likely, 
the flux of WR B is biased by being a part of a more extended 
emission region (A3; Sec. 3.3), a suggestion supported by data 
reduction that indicates the 8.6-GHz emission is resolved which 
is not expected for typical WR winds at this distance. For com- 
parison, the closest WR star to the Sun, WR 1 1 at ~0.26 kpc (van 
der Hucht et al. 1997) has a stellar wind that is resolved by the 
ATCA at 8.6 GHz, with a size of 0.47 arcseconds (Dougherty et 
al., in prep). If WR 1 1 was at a distance of 4 kpc, the wind would 
be ~ 30 mas in diameter and not resolved. This leads us to be- 
lieve that the source parameters derived for WR B are influenced 
by its location in the extended emission region A3. 

W44 (WR L) is the only WR star that has a spectral in- 
dex consistent with a partially optically-thick stellar wind, al- 
though the lower limit to the spectral indices of WR F leave 
this possibility open. For the WNVLh star WR L we adopt 
ju=2.0, Z=1.0 and 7=1.0 (e.g. Leitherer et al. 1997), resulting 
in M = 2 X 10-^(Voc/1000kms"^)/i/2 M^yi'^ respectively. For 
the WC9 star WR F adopting u=4.7, Z=l.l and y= 1.1 yields 
M = 3.3 X 10"5(v<«/1000kms"')/'^2 Moyr"'. Both these mass- 
loss rates are comparable to those found for similar stars in the 
Galactic Centre cluster by Martins et al. (2007) and other field 
WRs (e.g. Cappa et al. 2004). 

The remaining three WR detections (WR A, F, V) have flat, 
though highly uncertain, spectral indices. If they are flat, an in- 
terpretation due to composite spectra composed of both ther- 
mal and non-thermal emission components would be consistent 
with that observed in a number of WR+OB binaries (Dougherty 
& Williams 2000). Certainly, among WR stars with flat radio 
spectra, those that have been observed at high resolution (e.g. 



Williams et al. 1997; Dougherty & WiUiams 2000; Dougherty 
et al. 2005, and references therein), all have identified WCRs 
arising in CWBs. 

Recendy, Wd 1 has been observed by CHANDRA (Clark 
et al. 2007; Skinner et al. 2006; Muno et al. 2006a,b) and die WR 
stars WR A, B, F and L are aU X-ray bright (L^ ~ 10^2-33 erg s"^) 
and hard (kT > 2.6 keV i.e. T > 3 x lO'' K). These temperatures 
are expected for plasma heating by shocks in a WCR (Stevens 
et al. 1992; Pittard 2009; Pittard & Parkin 2009), which led both 
Clark et al. (2007) and Skinner et al. (2006) to suggest fliese WR 
systems are all CWB systems. In support of this argument, each 
of these WR stars have exhibited photometric variations in the 
optical (Bonanos 2007). WR B demonstrates clear eclipses while 
WR A demonstrates a 7.63-day periodic modulation to its light 
curve and is also found to be a spectroscopic binary (Crowther 
et al. 2006), though such a short-period system is not anticipated 
to exhibit evidence of non-thermal radio emission due to a wind 
collision (Dougherty & Williams 2000). However, if all these 
X-ray bright WR stars are interpreted as CWBs, then the appar- 
ently thermal radio spectrum of WR L would be attributed to 
a sufficiently short-period binary where any non-thermal radio 
emission from a WCR is completely absorbed by the circumbi- 
nary stellar wind(s) or the relativistic electrons are cooled by the 
intense UV radiation of a putative massive companion. 

The IR excess associated with WR F has been attributed to 
hot dust (Crowther et al. 2006; Groh et al. 2006), again sup- 
porting a CWB interpretation, since dusty WC stars are widely 
considered to be in binaries with OB companions where wind- 
wind interaction is the essential mechanism that attains the high 
gas densities required for dust formation (e.g. Williams 1999; 
TudiiU et al. 1999; Monnier et al. 1999; WilHams et al. 2009) 

A synthesis of the IR, X-ray and these radio observations 
therefore argues for a high binary fraction amongst the WR stars, 
with C08 suggesting that at least 70%, and possibly the complete 
population, of WR stars in Wd 1 are in massive binary systems. 
However, given the lack of radial velocity variations from spec- 
tral line observations this assertion has yet to be verified, and 
may be exceedingly chaUenging to verify. 



5. Extended emission 

Lastly, we consider briefly the extended emission and its ori- 
gin. Trivially, optically-thin thermal emission is associated with 
young massive clusters still embedded in Giant Hii regions, 
where the young massive stars ionize the remnants of the natal 
Giant Molecular Cloud. However, this seems unlikely for Wd 1. 
At an age of ~5Myr one might expect the action of a cluster 
wind to have sculpted a cavity ~ 10 parsec in extent e.g. the 
G305 star forming complex surrounding the young (>2-5Myr) 
massive clusters Danks 1 and Danks 2; Clark & Porter (2004), 
while the optically-thin thermal clumps associated with Wd 1 
are found significantly closer to, and indeed within the confines 
of, the cluster. Kothes & Dougherty (2007) report the presence 
of a large (25pc radius) H i bubble (their feature B3) in the neu- 
tral envirormient of Wd 1, which they attribute to the action of 
the cluster wind on the natal GMC and wider ISM. If the clus- 
ter wind swept the region clear of the natal material, the dif- 
fuse emission would then have to arise from material ejected 
recently from the steUar population of Wd 1, either via stellar 
winds or SNe. Assuming velocities of 20-200kms"' (appropri- 
ate for RSG/LBV and SNe ejecta) at a displacement of a ~parsec 
yields an age for material of 5xl03"^yr. Furthermore, a total ion- 
ized mass of 15 Mq only requires an average mass-loss rate of a 
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3 X 10"^ Moyr ' over the cluster lifetime, readily available from 
the population in Wd 1 as discussed in the previous sections. 

An origin in recent stellar ejecta seems to be certain for stars 
such as W20 and W26 with their large cometary nebula, but may 
not be so clear for other regions such as A6, where no associated 
stellar source is readily apparent in the FORS2 image. The ion- 
ized mass of A6 is estimated to be ~ 0.5 Mq, which could have 
certainly originated from a stellar wind. In addition, A6 does 
have a cometary-like structure that also appears to point towards 
the cluster centre region, like the nebulae around W20 and W26 
(see Fig. 4). An alternative explanation for some of the diffuse 
emission is that it arises in natal material left behind when the 
cluster region was largely swept clear by the cluster wind, pos- 
sibly through shielding by over-dense regions in the natal cloud. 
The nature and origins of the extended emission will be explored 
further in a forthcoming paper. 

6. Summary 

We have presented the results of a multi-frequency radio survey 
of the galactic Super Star Cluster Wd 1, and used optical, IR 
and X-ray observations from the literature to elucidate the na- 
ture of the radio sources detected. We detect 18 radio sources 
for which positional coincidence suggests an association with 
a cluster member. These radio stars are associated with every 
class of post-MS star present in the cluster. Moreover, they com- 
prise a diverse population of point-like, unresolved sources and 
extended, resolved sources, with spectral indices corresponding 
to thermal, non-thermal and composite thermal and non-thermal 
emission. Nevertheless, it appears possible to understand these 
properties under the current evolutionary paradigm for massive 
(binary) stars. In brief, the radio observations reveal the follow- 
ing properties for the massive stars in Wd 1 : 

- The brightest stellar radio source in Wd 1 is W9, with an 
8.6-GHz luminosity that places it amongst the most lumi- 
nous radio stars known. The emission comprises both a 
point source and a more extended nebula. The flux and 
spectral index of the point source implies an origin in a 
powerful stellar wind, with a mass-loss rate of 9.2^^5 x 
10-5(v„o/200kms"')/'/2 Mq yr"'. The extended nebula is 
deduced to arise from an earlier mass-loss epoch with a 
mass-loss rate ~3 times higher than the current stellar wind, 
close to the limit expected for a line-driven wind (Smith & 
Owocki 2006). An obvious comparison to make is to LBVs, 
with W9 having a current mass-loss rate similar to galactic 
examples e.g. the Pistol star. Likewise, the mass-loss rate de- 
duced for the extended nebula is comparable to several other 
galactic LBVs during outburst, although orders of magni- 
tude less than inferred for both P Cygni and rj Carinae dur- 
ing outburst (Clark et al. 2009, and references therein). The 
X-ray properties of W9 imply that it is a binary system, al- 
though it is not currently possible to constrain the nature of 
the companion. An alternative scenario with significant mass 
loss in a common envelope phase of stellar evolution cannot 
be excluded. 

- Surprisingly, three of the >100 evolved OB stars (luminosity 
classes Ill-Ia) are detected, with radio fluxes an order of mag- 
nitude larger than expected for stellar winds in these types. 
It is suggested they may be CWBs based on potentially com- 
posite spectra of both non-thermal and thermal emission. 
However, none of the three detected stars exhibit X-ray emis- 
sion characteristic of CWBs or evidence of RV variations due 
to a binary companion. Alternatively, their seemingly high 



radio luminosity may be due the influence of the extended 
emission in which they are embedded. 

- All four RSGs are detected, with three associated with large 
nebulae with ionized masses up to ~ 0.26 M©, emphasis- 
ing the importance of mass loss in this evolutionary phase. 
Of these, the nebulae around W20 and W26 have a pro- 
nounced cometary morphology, suggesting significant inter- 
action with either the intracluster medium or cluster wind. 
W237 shows less evidence for such interaction and has a 
kinematic age of ~3,600 yr and a time averaged mass-loss 
rate of 2 x 10^^(Voo/30kms"') Moyr"^ This is consistent 
with other field RSGs, although it is substantially lower than 
inferred for NML Cyg and VY CMa during the formation of 
their nebulae. 

- The YHG W4 is argued to have a stellar wind with a mass- 
loss rate of 10"^(Voo /200km s~^)f'^ Moyr"\ consistent with 
the few estimates available for other field YHGs. The ex- 
tended nebulae associated with W4, 12a and 265 are signif- 
icantly less massive than those associated with the RSGs in 
Wd 1. It is argued they arise from quiescent mass loss rather 
than during outburst episodes. 

- Since neither the YHGs nor RSGs are hot enough to ionize 
their own stellar winds and/or more extended nebulae, the 
requisite ionizing photons must arise from either an unseen 
companion or the cluster radiation field. 

- Of the extreme B-type hypergiants, only the LBV W243 was 
detected, with a spectral index consistent with thermal emis- 
sion. The corresponding mass-loss rate is directly compara- 
ble to that found for the YHG W4, as expected given the 
similarity in current spectral type and radio flux. Upper lim- 
its of 2 X 1 0-^(vo<, /200km s" ' Moyr"! for the three other 
B hypergiants were found, consistent with mass-loss rates 
amongst field stars of these types(Leitherer et al. 1995). 

- Five of the 24 WRs known in Wd 1 were detected. Of these, 
WR L has a partially optically-thick wind, with a mass-loss 
rate consistent with stars of identical spectral type in the 
Galactic Centre cluster and the general field population. The 
remaining three (WR A, B and V) are identified as having 
composite spectra from a CWB. The optical and X-ray prop- 
erties of WR A and WR B have previously indicated these to 
be binaries, while this is the first hint of binarity for WR V. 

In combination with X-ray, optical and near-IR datasets, 
these radio observations of Wd 1 provide striking evidence for 
the prevalence of binarity among the massive stars in Wd 1. This 
lends strong support to the previous estimates of high binary 
fractions, in excess of >40% (COS, Ritchie et al. 2009a). 

The rich stellar population of Wd 1 permits us to investigate 
the evolution of mass-loss rates as stars evolve from the MS. 
With the reduction of O-star mass-loss rales due to wind clump- 
ing (e.g. Repolust et al. 2004), it has been recognised that the ma- 
jority of mass loss must occur in the post-MS transitional phase. 
While the exact evolutionary sequence through the transitional 
'zoo' is uncertain for stars of M,„,7,v,/~40 M©, it is interesting that 
the radio mass-loss rates directly determined for the YHG W4, 
the LBV W243 and the WNLh star WR L, which form a direct 
evolutionary sequence in some schemes, are all closely equal to 
one another. 

However, with the exception of W9, such mass-loss rates 
(<lO"^M0yr~^) are likely insufficient to remove the H-rich man- 
tle unless stars remain in the transitional phase for significantly 
longer than expected. This in turn suggests that an additional 
mechanism is required to shed the requisite mass, with short- 
lived episodes of greatly enhanced mass loss an obvious candi- 
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date. Indeed mass-loss rates > 10 ^Moyr ' have already been 
inferred for RSGs (e.g. VY CMa; Smith et al. 2001) and directly 
observed for the YHG p Cas (Lobel et al. 2003). The nebulae 
around the RSGs W20, 26 and 236 already indicate that signif- 
icant mass loss has occurred for some stars within Wd 1 , while 
the mass-loss rate inferred for W9, over a magnitude greater than 
any other transitional star, is of obvious interest. Indeed, is W9 
undergoing an 'eruptive' event currently? 

We note that both the optical spectrum and spectral energy 
distribution of W9 indicate that it is not a cool hypergiant, sug- 
gesting that extreme mass-loss events are likely not confined to 
a single evolutionary state, such as RSGs, and instead occur in 
both hot and cool transitional phases. Indeed long term spectro- 
scopic observations indicate significant instability in the early- 
B hypergiantAVNVL, RSG and YHG populations (Clark et al. 
2010). Hence one may not a priori assume a single physical 
mechanism leads to enhanced mass loss in post-MS stars. When 
combined with detailed modeling of the stellar spectra, further 
observations to derive the nebular properties of cluster members, 
such as expansion velocity and chemical abundances, would be 
very valuable in order to constrain the phases in which enhanced 
mass loss occurs; thus Wd 1 will be invaluable for investigating 
the physics that mediates the passage of MS O-type stars to the 
Wolf-Rayet phase. 
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Appendix A: Modelling radio emission from a 
circumstellar envelope 

For a spherically symmetric isothermal plasma envelope of tem- 
perature T and outer radius Ro at distance D from the Sun, the 
total flux at frequency v is given by 



„-rv(P) 



)P dp , 



(A.l) 



where Ty(p) is the optical depth along a line-of-sight at distance 
p from the line-of-sight to the centre of the envelope^* (Fig. A. 1) 
For a Maxwellian plasma in LTE the opacity is given by 
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for electron temperature Tg, gv flie gaunt free-free factor (see e.g. 
Leitherer & Robert 1991), Z the rms charge of the plasma , y the 
ratio of electron to ions, and n,, is the electron density along a 
line-of-sight segment length dz (Rybicki & Lightman 1986). 

For a steady-state, constant velocity stellar wind of mass- 
loss rate M and terminal velocity Vco, the equation of continuity 
const, leading to 
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with fi the atomic weight of the plasma, and m// the mass of 
hydrogen. The free-free opacity is then given by 
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Analytic solutions for the integral hip) in Eqn. A.2 and the en- 
velope geometry shown in Fig. A. 1 can be derived, namely: 
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Fig. A.l. Geometry of a model circumstellar envelope with outer radius 
Rg and inner radius Rj. For Rj < r < R„ the plasma density oc r^^ and is 
zero for r < Rj. 



if Ri < p < Ro- 

Determining the flux from an extended nebula at any given 
frequency requires using Eqns. A.2 and the two conditional cases 
for h)p) given by Eqn. A. 3 and A.4 to determine the opac- 
ity along a particular line-of-sight, before solving Eqn. A.l. To 
model the flux from a "standard" stellar wind, R, is set to the 
radius of the underlying star Rt (for the models considered here 
Rt « Ro, so it is assumed that = 0) and only the conditional 
case for hip) given in Eqn. A.4 needs to be considered. 



' p is often referred to as the "impact parameter". 



